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Abstract. We study in details a pumping mechanism for the A = 1.35 cm maser transition 6i6 — > 523 in ortho-KkO 
based on the difference between gas and dust temperatures. The upper maser level is populated radiatively through 
4i4 — * 5os and 5os — > 6i6 transitions. The heat sink is realized by absorbing the 45 /im photons, corresponding to 
the 523 — * 4i4 transition, by cold dust. We compute the inversion of maser level populations in the optically thick 
medium as a function of the hydrogen concentration, the gas-to-dust mass ratio, and the difference between the 
gas and the dust temperatures. The main results of numerical simulations are interpreted in terms of a simplified 
four-level model. We show that the maser strength depends mostly on the product of hydrogen concentration 
and the dust-to-water mass ratio but not on the size distribution of the dust particles or their type. We also 
suggest approximate formulae that describe accurately the inversion and can be used for fast calculations of 
the maser luminosity. Depending on the gas temperature, the maximum maser luminosity is reached when the 
water concentration Nh 2 o ~ 10 6 -j- 10 7 cm -3 times the dust-to-hydrogen mass ratio, and the inversion completely 
disappears at density just an order of magnitude larger. For the dust temperature of 130 K, the 6ie — > 523 
transition becomes inverted already at the temperature difference of AT ~ 1 K, while other possible masing 
transitions require a larger AT > 30 K. We identify the region of the parameter space where other ortho- and 
para-water masing transitions can appear. 
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1. Introduction 

Sources of strong water maser emission at wavelength 
A = 1.35 cm have been found in many astrophysical ob- 
jects such as active galactic nuclei, carbon rich stars, pro- 
tostellar regions and comets. Maser emission is a powerful 
tool for investigating the physical conditions in the emit- 
ting regions because of its high brightness as well as a 
high sensitivity to the physical parameters of the medium 
in which the maser amplification takes place. 

There are several mechanisms for the ortho-H^O 6i6 — > 
523 (A = 1.35 cm) maser pumping, de Jong (1973) pro- 
posed a comprehensive model where the breakdown of 
thermal equilibrium occurs in the surface layers of the 
dense gas cloud. Upon approach to the surface, the 523 ~~ ► 
4i4 transition (45 fim) (which is among the most impor- 
tant for the maser action, see Fig. 1) becomes transpar- 
ent first and level 523 becomes underpopulated. The de 
Jong model was widely used to describe the water masers 
from late-type stars and star-forming regions (Cooke & 
Elitzur 1985; Elitzur, Hollenbach, & McKee 1989; Neufeld 
& Melnick 1991), and the molecular accretion disks in ac- 
tive galactic nuclei (Neufeld, Maloney, & Conger 1994; 
Babkovskaya & Varshalovich 2000). 



The de Jong mechanism is, however, not able to ex- 
plain the most powerful masers in star-forming regions. A 
different physics has to be involved. In shocks, the depar- 
tures from the equilibrium could be produced by collisions 
with hotter "superthermal" hydrogen gas (Varshalovich 
et al. 1983). Strelnitskij (1980, 1984) proposed that colli- 
sions with two species (e.g. charged and neutral particles) 
of different temperatures and with comparable collision 
rate can produce the necessary inversion at high hydro- 
gen densities needed to explain the high observed maser 
luminosities. Kylafis & Norman (1987) pointed out that 
the magnetohydrodynamic shocks produce naturally con- 
ditions where the electron temperature T is larger than 
the hydrogen temperature Tjj. In such a case, the lower 
maser level becomes underpopulated, because the rela- 
tive importance of collisions with neutrals is larger for the 
transitions to/from this level. The inversion in this model 
appeared only in a narrow range of the ionization frac- 
tion 10~ 5 - 10~ 4 and for T H < 60 K. At high ionization, 
electrons start to dominate the collisions and the levels 
are thermalized at T e , while at lower ionization, neutrals 
are dominating and the thermalization happens at Th- 
Using better estimates for the collision rates, Elitzur & 
Fuqua (1989) and Anderson & Watson (1990) ruled out 
this model unless the neutral particles are hotter than the 
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charged ones, conditions that are difficult to imagine in 
any astrophysical environment. 

All astrophysical objects that show water maser emis- 
sion are also expected to have non-negligible quantities 
of dust. If the dust and the gas have different temper- 
atures, the departures from the equilibrium are possible 
(Goldreich & Kwan 1974; Kegel 1975; Strelnitskij 1977; 
Bolgova et al. 1977). Goldreich & Kwan (1974) suggested 
that the radiation from the hot dust excites the water 
molecules to the vibrational state, and the heat sink is 
provided by collisions with cooler (than dust) hydrogen 
molecules. The possibility of the inversion in such a situ- 
ation was questioned by Deguchi (1981) who pointed out 
that because the collisional de-excitation rate between vi- 
brational states is much smaller than the pure-rotational 
collisional rate, the rotational collisional thermalization 
may quench the maser when vibrational collisional de- 
excitation becomes dominant (see also Strelnitskij 1988). 

Alternatively, the cold dust can produce the necessary 
inversion (Strelnitskij 1977; Bolgova et al. 1977). In the 
optically thick environment, the excitation temperature 
takes the values between the dust and the gas tempera- 
tures depending on the relative role of dust and collisions 
in the destruction of the line photons. Deguchi (1981) con- 
sidered the following cycle of the maser levels pumping 
4i4 — ► 5o5 — > 6i6 — ► 523 - * 4i4 (see Fig. 1). He showed that 
the downward transition 523 ~~ * a t 45 /xm is much more 
affected by the dust absorption than the upward transi- 
tions at A <~ 80 — 100 tun, because there is a strong peak 
near 45 /zm in the absorption coefficient of the cosmic- 
type ices (e.g. Moore & Hudson 1994). The upper level 
excitation temperature is then close to the gas tempera- 
ture, while the lower level becomes populated at the dust 
temperature. Similar inversion occurs for other types of 
dust too. 

One should note that in this model an arbitrary thick 
layer can participate in the maser action provided the gas 
and dust temperatures sufficiently differ. Chandra et al. 
(1984a) computed the maser efficiency when both the sur- 
face escape and the cold dust absorption mechanisms op- 
erate together. Recently Collison & Watson (1995) and 
Wallin & Watson (1997) applied this model to the masing 
disk in the active galaxy NGC 4258 concluding that the 
cold dust model is much more efficient (see also Neufcld 
2000). Thus the cold dust-hot gas model seems to be the 
most promising to explain powerful water masers in many 
astrophysical objects. The temperature difference between 
the gas and the dust can appear as a result of shock heat- 
ing or illumination by the UV- or X-ray photons, and/or 
because of the presence of the dust particles of different 
types and sizes which assure different temperatures. 

The purpose of the present study is to determine the 
maser strength for a broad range of the physical parame- 
ters: the hydrogen concentration, the water-to-dust mass 
ratio, the gas and the dust temperatures as well as the 
dust type. Previous studies (Deguchi 1981; Chandra et al. 
1984a) used the rectangular line profile instead of the 
Doppler one as well as the collisional coefficients with 
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Fig. 1. A portion of the diagram of the ortho-HbO showing 
rotational levels involved in the pumping of 6i6 — ► 523 maser. 



hydrogen of Green (1980), which later have been much 
improved (Green, Maluendes, & McLean 1993). Only a 
handful set of parameters was explored in the more re- 
cent investigation of Collison & Watson (1995) and Wallin 
& Watson (1997) who considered the case of saturated 
maser. 

In its full statement this problem needs the simultane- 
ous solution of the statistical balance equations together 
with the radiative transfer equations for all spectral lines. 
However, because the cold dust - hot gas model can oper- 
ate inside a molecular cloud where almost all transitions 
(except masing) are optically thick, the radiative transfer 
can be handled in a very simple manner using the escape 
probability formalism, where now the role of the escape of 
spectral line photons is played by the dust absorption. 

The paper is constructed as follows. We give the formu- 
lation of the problem and describe the numerical method 
in Sect. 2. We develop a simple four- level model that con- 
tains most of the physics involved in Sect. 3. The numer- 
ical results for different dust types and size distributions 
are presented in Sect. 4, where we also propose simple 
analytical formulae for the inversion which describe the 
results of simulations. Finally, in Sect. 5 we present our 
conclusions. 



2. Method 

2.1. Radiative transfer 

We consider a slab consisting of a mixture of molecular 
hydrogen, water molecules and dust particles. The radia- 
tive transfer equation for the specific intensity I x describ- 
ing the transfer of the line radiation in the presence of 
the continuum absorption and emission is (Hummer 1968; 
Mihalas 1978) 

A*-r^ = -{u-L<t> x + (Xd)I x + a L 4>xS u i + a d B u i(T d ), (1) 
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where z is the geometrical depth, fi is the cosine of the an- 
gle between the direction of propagation and the outward 
normal, 



c | n l 9u , 
Oul = 1 

n u gi 



, u > I, 



(2) 



is the source function for the line transition between the 
energy levels u and I, and B u i(T d ) is the Planck function 
at the transition frequency v u i characterized by the dust 
temperature. All the intensities and the source functions 
are given in units 2hv^ ll /c 2 . The dust absorption coefficient 
ay = K dPd depends on the frequency-dependent (constant 
within the line) opacity n d [cm 2 g _1 ] and the dust density 
Pd = /d2m p A^H 2 [g cm -3 ], where m p is the proton mass, 
_/Vh 2 is the concentration of molecular hydrogen, and f d is 
the dust-to-hydrogen mass ratio. 

The line averaged absorption coefficient (includ- 
ing induced emission) depends on the corresponding level 
populations: 



«l = g u A u iAn ul 
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87Tc(Af D /f M; ) 
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(3) 



where An u i = ^ — 2a., n u is the fractional level pop- 
ulation (J2 u nu = 1)' ^h 2 o is the water concentration, 
g u is the statistical weight for a given level, A u i is the 
Einstein coefficient, and A„; = c/v u i is the transition wave- 
length. We further assume complete frequency redistribu- 
tion within the line, i.e. both the absorption and emission 
line profile are described by the same function <f> x normal- 
ized as J 4> x dx = 1, where x = (v — u u i)/Avn is the fre- 
quency within the spectral line in thermal Doppler units 
A^d = v u i(2kT /mc 2 ) 1 / 2 . For the Doppler profile assumed 
here, <p x = tt^ 1 / 2 cxp(— x 2 ). 

Let us note the averaged over line profile optical thick- 
ness of the slab in u/-line as 2tq = J a^dz. A formal 
solution of the radiative transfer equation (1) for a homo- 
geneous (i.e. T — const, T d = const) slab gives the mean 
intensity of the radiation averaged over the line profile at 
the optical depth r (Hummer 1968; Mihalas 1978): 

//■2ro 
<j> x J x (r)dx = (1-6) J Sm^K^t - r)dt 

/•2to 

+ 5B ul (T d ) / Li(t-T)dt, (4) 
Jo 



where the kernels 
1 



KM 



E.a^+pMcpidx, 



2(1-6) 

Li(s) = E 1 ([cf> x + (3}\ S \) ( t> x dx 



(5) 
(6) 



are normalized as Ki(t)dt — Li(t)dt = 1. Here 

E n (r) — f u n ~ 2 exp(— r/u)dfi is the exponential integral 
function, 



(3 = «d/aL, 



(7) 



and 



6 = (3 



b x + (3 



dx 



(8) 



(index ul is omitted in (3 and 6 as well as in the formulae 
below) is the probability per single interaction act that 
the line photon will be absorbed by dust (obviously, 1 — 5 
is the probability for a photon to excite a molecule). 

If the source function does not vary much, we can take 
it out from the integral (the so called on-the-spot or the 
first order escape probability approximation) and get: 



J(r) ~ (1-5) 



1 - \k 2 (t) - \k 2 (2t 
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r) 



r) 
B(T d ) 



S(r) 



(9) 



(l-p)S+p c B(T d ) 



where 



/oc- 
JTi(t)dt 

= ihJ^h E ^ + ^ dx > (10) 

/oc 
Li(t)dt 

For every column density of the slab we can take the inten- 
sity in the middle plane J(to) as the representative one. 
We perform all our calculations using equation (9) with 
r = r = a^H (where H is the slab half-thickness), i.c 
taking 1 



P = 5+(1-6)K 2 (t), P c = 6[1 - L 2 (t)}. 



(12) 



If, however, the line is optically thick, the escape is negli- 
gible and the mean intensity is simply 



J~ (1-6)S + 5B(T d ). 



(13) 



This expression can be obtained directly from equation 
(1) assuming isotropic intensity of radiation (i.e. assuming 
dl x /dz = and I x = J x ), solving for J x , and integrating 
it over the frequency with a weight <j) x . If most of the con- 
sidered lines are optically thick, the populations do not 
depend on the optical depth of the considered point but 
are determined by the local conditions (such as hydrogen 
density, gas-to-dust density ratio, and dust and gas tem- 
peratures) only. 

Following Hollenbach & McKee (1979), we approxi- 
mate 5((3) as 



2(3 



1 + 2(3 



In e + 



1 



^(3 



1/2 



(14) 



1 Collison & Watson (1995) and Wallin & Watson (1997) 
have considered a semi-infinite slab with To = oo and used the 
escape terms p = 8 + (1 — 5)K2(t) and p c = p — K^((3 = 0, t), 
while the correct ones in that case are p = 8 + (1 — 8)K% (t)/2 
andp c = S[l -L 2 (t)/2]. 
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and slightly modify the expression for K2(t): 
cxp(-T d ) 1 



K 2 (r)~ 



1 + (1 - 5)V^T l + r c [2^1n(2.13 + r c 2 )] 1 /2 



(15) 



where r c = r/y/ir is the optical depth in the line center 
and Td = (3t is the optical depth for the dust absorption. 
The continuum escape coefficient p c = 5[1 — L 2 {t)] is ap- 
proximated as 



p c ~ X [5 - 



exp(-r d ) 



l+T c [TT\n{1.2b + T c )Y/* 

S + T d [l-ln(l + l/T,)Y 



1+T d 



(16) 



X = 1-0.095- 



1 + y 



•ln(l + 0.028//3), 



where y = (0.477t c ) - 277 , t* = max(rd,T c ). These expres- 
sions for K 2 and p c are accurate within 10% for any (3 
and r. 

The masing lines need a different treatment. We ne- 
glect the dust influence in such transitions, i.e. assume 
ft = 5 = p c = 0. Using expression (10), one can show that 
at large maser optical depth |t| 3> 1 



exp(|r c |) 



:|\A"ln| 



(17) 



Then the line escape coefficient for r < can be repre- 
sented as 



p(r) 



exp(|T c |) + 1.18|r c | 
1 + |tc| [ttMI + ItcI)] 1 ^- 



(18) 



This expression reproduces the exponential asymptotic 
behavior at large |r c | (17) as well as matches the asymp- 
totic of K 2 at small r c (see eq. 15). 

We consider different types of dust. The amorphous ice 
absorption coefficient is computed from the optical con- 
stants of Leger et al. (1983). The crystalline ice absorp- 
tion coefficient are based on the data from Bertie et al. 
(1969) and Warren (1984). The silicate and graphite data 
are from Laor & Draine (1993). 

2.2. Statistical balance equations 

The populations of M levels can be determined from 
the population balance equations which, in the station- 
ary case, take the form: 

^2mW lu = n u ^W uU u=l,2...,M-l, (19) 

M 

^2 n u = i, 

where W u i = C u i + R u i is the total rate of the u — > I 
transition, 



rA 9u^i ( hv ul 



Ci = N U2 k ul (T) (20) 



being the rates of the collisional excitation and deexci- 
tation of water molecules by molecular hydrogen which 
kinetic temperature equals that of the water. The rates of 
the radiative excitation and deexcitation are given by 



R l ul = A ul (l + J u i) . 



(21) 



The Einstein A-coefficients are taken from Chandra et al. 
(1984b). Collisional deexcitation rates are from Green, 
Maluendes, & McLean (1993), while the excitation rates 
are computed using the detailed balance condition. We 
take into account the first 45 rotational levels of ortho- 
and para-H20 molecules and all possible radiative and 
collisional transitions between them. 

Substituting the radiation intensity from (9) into (21) 
and the statistical balance equations (19), we obtain: 



^ A lu (p lu ni + pf u Bf u giAn u i) 

1>U 

^2 A u i (Puin u + p C uiB^g u Aniu) 



(22) 



Ku 



= ^2 (Cuiriu - Ci u ni) , u = l,2,...,M-l, 

with the same normalization as before ^2 U n u = 1. 

The populations can be found from the solution of 
the system of non-linear algebraic equations (12), (18), 
and (22). We use the Newton- Raphson method with line 
searches and backtracking (Press et al. 1992). We nor- 
mally start from high hydrogen density where the solu- 
tion is close to the Boltzmann at the gas temperature and 
use the obtained solution as a zeroth approximation for 
the next point. We continue iterations until the maximum 
error in the system becomes smaller than lO^ 11 and the 
maximal change in populations is smaller than 10~ 9 . It 
takes 5 to 15 iterations for most sets of parameters and 
each computation takes on average about 15 ms CPU time 
on a Pentium IV 2 GHz Linux PC. 

2.3. Maser intensity 

If the maser amplification takes in the homogeneous 
medium and there is no velocity gradients, the maser in- 
tensity can be estimated as 



4(Vm) - (h + l^mDexp^m^ 



(23) 



where Iq is the background continuum intensity, S m is the 
source function for the u — > I masing transition, and the 
optical depth in the maser line T m is 



^m-^coh — 



x li9uA u i 
8irc(Avr> / Vui) 



An m Nn 2 oL co h, 



(24) 



where An m = n u /g u — ni/gi (positive for the maser) 
and L co h is the typical coherent length of the maser am- 
plification. If the excitation temperature of the masing 
transition is larger than the brightness temperature of 
the background, one can neglect Io (and obviously other 
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way around too). We assumed the first possibility (i.e. 
Jo =0). In this paper we quote values of r m for L co h = H. 
Obviously at grazing angles to the slab surface this value 
can be exceeded many times. 

2.4. Parameter range 

The water level populations depend on the hydrogen den- 
sity _/Vh 2 , the dust and water concentrations and tempera- 
tures as well as the slab half-thickness H. When the pho- 
ton escape from the surface is negligible because the slab 
is optically thick either in continuum (jd > 1) or in the 
line (t> 1), the terms K 2 and L 2 can be omitted in the 
escape probabilities (12) (i.e. p = p c = 5). The solution 
of the population balance equations then depends on 5, 
which in turn depends only on the ratio of the absorp- 
tion coefficients (5, but not on the absorption coefficients 
individually. Since the dust opacity k& has a rather weak 
dependence on the size of the dust particles, it is rather 
natural to use the water-to-dust mass ratio /h 2 o//c1 as 
a parameter, where /h 2 o is the water-to-hydrogen mass 
ratio and /a is the dust-to-hydrogen mass ratio. The wa- 
ter concentraion relates then to these parameters through 
Nn 2 o = -^H 2 (/H 2 o//d)/d/9 (factor 9 comes from the ra- 
tio mn 2 o / mn 2 ) . When the escape through the surface is 
small and masers are unsaturated (i.e. a masing transition 
does not affect the level populations), the populations do 
not depend on /d and the scale-height H, and the opti- 
cal depth is then linearly proportional to fdH. In such a 
situation the main parameters are Nn 2 , /H 2 o//d, T and 
T<j. If the dust is suffiently cold, its own radiation is then 
negligible and dependence on Td disappears. We normally 
consider the dust colder than the gas, AT = T — Td > 0, 
because this is required for inversion to occur (in the op- 
tically thick case). 

One can also point out that for large fn 2 o/fd when 
(3 ~ 5 <C 1 in the main infrared pumping lines, the solu- 
tion of the balance equations (22) depends on the product 
^H 2 /H 2 o//d, but not on N H2 and / H2 o//d individually, 
because C u i/p u i oc Nn 2 /f3 oc iV H2 /H 2 o//d- 

We consider a range of the hydrogen concentrations 
from 10 8 to 10 125 cm -3 . At higher Nn 2 , the inversion is 
absent because of the thcrmalization by collisions with the 
hydrogen molecules, while at smaller Nn 2 we can extrap- 
olate the results from JVh 2 = 10 8 . 

In astrophysical objects it is easier to estimate the 
dust-to-gas mass ratio /d than other parameters of in- 
terest. We fix it at a standard for the interstellar medium 
value of /d = 0.01. This then determines the dust optical 
depth 



Td = 0.033«d 



^Vh 2 /d H 
10 10 0.01 10 14 ' 



(25) 



For the majority of calculations we consider the dust 
in the form of the amorphous ice grains of the size a = 
0.1 /im. Comparison is also made to different sizes and 
temperatures of dust as well as other types of dust (crys- 
talline ice, silicate, and graphite). 
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Fig. 2. Dependence of the 6i6-523 maser absorption coefficient 
im/(10 _14 cm _1 ) on the half-thickness of the slab for six sets 
of parameters (Nh 2 , fa 2 o/fd): (1) 10 9 ,1(T 2 ; (2) 10 10 ,10~ 3 ; (3) 
10 n ,10" 4 ; (4) 10 9 ,1(T 4 ; (5) 10 10 ,10~ 5 ; (6) 10 n ,10~ 6 . Other 
parameters are fixed at T d = 130 K, T = 250 K, f d = 0.01. At a 
slab thickness less that 10 cm, the effect of the boundary be- 
comes significant. For higher H, the maser becomes saturated 
when r m = 30. 

The water-to-dust mass ratio probably varies by or- 
ders of magnitude from one object to another. Therefore, 
we consider a broad range of fn 2 o/fd from 10~ 8 to 10°. 
The maximum possible /h 2 o m the interstellar medium is 
about 10 -3 which transforms to the maximum /h 2 o/ fd ~ 
0.1 for /d = 0.01. However, for a smaller dust content, 
/H 2 o//d can be larger. 

Figure 2 shows the absolute value of the absorption 
coeffient (in units of 10 _14 cm _1 ) in the main masing 
line 616 — ► 523 f° r six representative sets of parameters 
(Nn 2 , /H 2 o//d)- One sees that at H < 10 13 cm, the sur- 
face escape of photons (the de Jong mechanism) starts to 
influence the populations. At high H , the optical depth r m 
can be sufficiently large to saturate the maser. One can 
easily estimate the optical depth when it happens. The 
saturation occurs when the induced transition rate A J in 
the masing line becomes comparable to the rates of tran- 
sitions from the upper masing level in other strong lines 
(see eq. 21) which are R u i ~ 1. Because the Einstein A 
coefficient in the 1.35 cm 616 — > 523 masing line is about 
10~ 9 (and about 10~ 5 in other masing transitions), the 
saturation occurs when J <~ 10 9 . On the other hand, since 



J ~p\S\~exp(T m /V^)\S\ 



(26) 



(see eq. 18) and |5| <~ 30, the maser saturates at the op- 
tical depth r m ~ 30 (which weakly depends on the value 
of S). This limiting optical depth is reached at different 
H depending on the conditions. Here the maser saturates 
and the inversion decreases. 

The inversion is rather flat around H <~ 10 14 cm, where 
it is a function of the local conditions only. Therefore, we 
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Fig. 3. Levels of the constant maser absorption coefficient 
am/(10 _14 cm _1 ) (solid curves) and the maser optical depth 
r m (dotted curves) on the plane slab half-thickness - gas tem- 
perature for 7Vh 2 = 10 9 cm -3 and /H 2 o//d = 10~ 2 . Other pa- 
rameters are fixed at T d = 130 K, / d = 0.01. 



use this height for our calculations. We should also note 
that our results do not depend on the geometry of the 
system, and can be applied not only to the slab but to 
any other geometry. 

The contours of the constant maser absorption coeffi- 
cient and the optical depth on a plane H - gas temperature 
T are shown in Fig. 3. One sees that, independently of the 
gas temperature, the absorption coefficient is rather flat 
around H ~ 10 14 cm. The flat part becomes shorter at 
higher gas temperature, because the saturation (r m ~ 30) 
happens for smaller H. One should note here that the pa- 
rameter set ./Vh 2 , /h 2 o I Id taken at this graph gives much 
shorter flat parts than other sets presented in Fig. 2. Our 
fiducial H = 10 14 cm still seems to be a good choice for 
studying in details the physics of the unsaturated water 
maser in a dusty environment. It is also interesting to note 
that for small H inversion exist even when T < Td because 
of the action of the de Jong mechanism. 



3. Analytical four level model 

Our main goal is to study in details the inversion mecha- 
nism of the 6i6 — > 523 masing transition. Before we pro- 
ceed to the results of calculations, it is useful to under- 
stand the physical processes responsible for the action of 
this main maser using a simplified system of molecular 
levels, which mainly participate in the maser pumping: 
4 14 (level 1), 5 05 (level 2), 5 23 (level 3) and 6 16 (level 4) 
(see Fig. 1 and Deguchi 1981). The upper maser level 4 
interacts mostly with the 2nd level which in turn interacts 
with level 1, while the lower maser level 3 interacts directly 
mostly with level 1. (We neglect here rather strong tran- 
sitions from the 2nd to the 3rd level via level 5i4.) The 



rates of radiative (in the case of the unsaturated maser) as 
well as collisional transitions via the masing line 4 — > 3 are 
much smaller than the rate of transitions to other levels. 
Thus, we can assume that the populations at the masing 
levels are completely unrelated to each other, and depend 
only on the rate of the exchange to other levels. This al- 
lows us to write a system of three (two-level) population 
balance equations that relate the corresponding popula- 
tions in the standard form: 



n u {R u i + C u i) = ni(Ri u +C iu ), 



(27) 



where ul are 21, 42 and 31. From equation (2), we get the 
corresponding (dimensionless) source functions: 



S = 



1 7 + 



1 + e 



1 + 6 



-B, 



(28) 



where e = C(l — exp[— E/T])/A, E is the corresponding 
transition energy in Kelvin and we omitted the indices ul. 
Assuming that the photon escape is negligible, we substi- 
tute the expression for the radiation field (13) and obtain: 



-B 



1 



i + e i+£ 



-B a 



where £ = e/5. The intensity is then 



-B 



1 + 



i+e i+e 



(29) 



(30) 



This is, of course, not a self-consistent solution, because S 
depends on j3 which in its turn is a function of the level 
populations. We can assume in the first approximation 
(for calculating (3) that there is a Boltzmann distribution 
of the populations at the gas temperature T. We thus see 
that the source function is determined by a single param- 
eter, the ratio £. We also can note that for a large water 
content (or small dust content, i.e. /3<g;l)<5^/3ocl/^4 
and £ does not depend on the line strength (since A u \ can- 
cels out). Then S is completely determined by the relative 
ratio of the photon destruction rate by collisions to that by 
dust absorption. On the other hand, when the dust is dom- 
inating the absorption (i.e. /3 » 1), 5 <~ 1 and £ = e, the 
source function is determined just by the relative impor- 
tance of collisional and spontaneous radiative transitions. 
Now we can obtain the general expression for the inver- 
sion in the levels 4-3. Since n u /g u = ni/gi S u i/(l + S u i), 
the inversion is 



An m = 



ni 
9i 



S 2 



S42 



'31 



1 + Sal 1 + S42 1 + S 31 



(31) 



Because the 1 <-> 3 transition is much more affected by 
dust that other main transitions, it is possible that the 
population of the lower masing level is determined by the 
dust temperature (i.e. S31 = -B 31 ), while in other transi- 
tions the dust influence is still negligible (i.e. S42 = B42, 
and 52i = B 2 i). The inversion then reaches the maximum 



9i v 



-E41/T 



-E 31 /T c 



)■ 



(32) 



Let us now consider some limiting cases. 
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3.1. Large water content 

When the water concentration is large, the source func- 
tions depend on one parameter £ = e/S. This means that 
the inversion is the function of A^/i^o/Zd only. If £ for 
all main lines is very small (e.g. small hydrogen density) 
then the levels are thermalized at the dust temperature, 
while for large £, thermalization occurs at the gas temper- 
ature. The limits on Ah 2 where the inversion is possible 
depend on fa 2 o/fd- 

Let us first investigate the limit ( « 1. The source 
functions can be represented then as S = B d +^AB, where 
AB = B — B d . Then 



1 + 5 



= cx P (-£yr d ) 



AB 



B d (l + B d ) 



(33) 



For the amorphous ice the £ parameters can be approxi- 
mated as 



&l « 6 x 10- 5 iVH 2 /H 2 o//d, 
£42 w 10~ Nu 2 fu 2 o/ fd, 

&1 « 10- 7 AT H2 /H 2 o//d- 



(34) 



Because, £31 <C £42 <C £21 > we can keep only the terms with 
£21 in equation (31). A condition of the positive inversion 
An m > then transforms into 



£21 > £21, min — 



E m B d 1 (l + B d 1) 



21 



(35) 



which corresponds to iVH 2 /H 2 o//d > 120 for T<j = 130 
K. Here E m = E41 — E31 w 1 K is the maser transition 
energy. For smaller amount of water, the inversion disap- 
pears. However, we implicitly assumed here that £AB is a 
small correction to B d . If the dust is sufficiently cold, this 
is not true anymore. Thus, if B d <C £B <C 1, the source 
function is just S w £B. Because £21 is by far the largest 
among the considered transitions, the inversion appears 
when 



61 > £21, 



31 



31 



£42 B21B42 ' 



(36) 



corresponding to Nn 2 fn 2 o/fd > 50. Thus generally the 
inversion exists at £>i > max[£ 2 i,min, Sn.miJ- 

In another extreme case, small influence of the dust 
relative to collisions, i.e. £ » 1, we get S = B — A£?/£. 
Now keeping only the term oc 1/^31 we obtain a condition 
for the inversion 



^31 < £31, max 



AH 



31 



E m £31(1+531)' 



(37) 



which corresponds to Ah 2 /h 2 o/ fd < 10 9 for T = 250 K. 

We can get an estimate for the temperature difference 
needed to produce the inversion. The inversion appears 
first in the region where the 3 — > 1 transition is dominated 
by dust (i.e. £31 < 1 and S = B d + £AB), while for other 
main transitions the dust influence is still negligible (i.e. 
£21 > £42 > 1 and S = B - AB/£). Expanding AB in 



the vicinity of T d as AB = (B d /T d ) 2 EAT, wc get from 
the condition An m > 0, 

E m T<\ 



AT > 



E31 (1 - fcie-^i/Td) - E 42 e- E ^/ T */U 2 ' 



(38) 



where we neglected the terms of the order 1/^21 • For 
Nn 2 fn 2 o/fd = 10 6 , we have £31 = 0.1 and £ 42 = 10, and 
then AT > 0.5 K for T d = 130 K. In the limit £31 < 1 
and £42 S> 1, corresponding to S42 = B42, S21 = B21 and 
S31 = £?3i, we get a much simpler expression 



AT > E m T d /E 31 = 0.003T d . 



(39) 



We see that already a small difference in the temperatures 
produces the inversion. 

3.2. Small water content 

A small water content (or large dust content) is described 
by [3 > 1 and 6 ~ 1. For the slab thickness H = 10 14 
cm, the medium is transparent for small fn 2 o/fd and 
small Ah 2 (since the dust optical depth is proportional 
to the hydrogen density), and thus a similar analytical 
description as above is not valid. However, for larger H, 
and/or Nn 2 , and/or / d , the medium still can be suffi- 
ciently opaque. Then, the source functions in all pumping 
lines are 

; (40) 



S 



1 



-B + 



1 



-B d . 



Thus the inversion does not depend anymore on /H 2 o//d, 
but is defined only by the hydrogen density (which deter- 
mines e) and the temperatures. At very high Ah 2 ,£>1 
and all source functions are given by the Planck function 
at the gas temperature, while for small Ah 2 , e < 1 and 
thermalization occurs at the dust temperature. We can 
determine the condition when the inversion appears from 
the equation An m = 0. For small Ah 2 , all e <C 1 and 
the source functions for every transition can be written 
as S = B d + eAB. Substituting this into equation (31) 
and keeping only the terms of the first order in e, one gets 
the lower limit on the density where the inversion is still 
possible: 



^min = E™ [ AB 2 1 (e 2 l /Ah 2 ) 

H2 " T d 1 B^l + Bi,) 



(41) 



+ 



AB 42 (642/Ah 2 ) _ Ag 3 i(e 3 i/AH 2 ) 
B 4 d 2 (l + B 4 d 2 ) B 3 d i(l + £ 3 d i) 



(note that e is proportional to the hydrogen density). For 
example, for T = 250 K and T d — 130 K, we get limiting 
Af 2 in = 2 x 10 8 cnT 3 . 

Analogously, we can get an upper limit on Ah 2 assum- 
ing that e » 1. Rewriting S = B — AB/e and keeping the 
terms of the first order in 1/e we get 



r ax = T_ f Aj? 31 (A H2 /e 3 i) 
H2 ' E m \ B 31 (l + B 31 ) 

AB 2 i(A H2 /e 2 i) AB 42 (A H2 /e 42 ) 



(42) 



-821(1 + B21) -842(1 + 542) 
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Fig. 4. Contour plots of the levels of the constant inversion An m (dotted curves) and constant optical depth (solid curves) 
for the 6i6~523 transition at the plane hydrogen concentration - gas-to-dust concentrations ratio. A slab of the half-thickness 
H — 10 14 cm and the dust in the form of the amorphous ice of the size a — 0.1/im was assumed with the dust-to-gas mass ratio 
fd = 0.01. The panels correspond to the following temperatures (a) T d = 30 K, T = 250 K; (b) T d = 130 K, T = 132 K; (c) 
T d = 130 K, T = 250 K; (d) T d = 130 K, T = 400 K. The dashed contour in panel (b) corresponds to the negligible photon 
escape case p — p c — 6. 



For the same parameters as above, this expression gives For our fiducial Td 



iV™ ax - 10 13 cm- 3 . 

Now let us obtain a minimum temperature difference 
needed for the inversion. Representing the source function 

as 5" = B d + -r^—AB and expanding AB for small AT, 



we get 



AT, 



1 + e 



Em T d < fHi^i c -£ 21 /T d 



1 + C21 



£ 42 T 42 Ei2 / Td _ £31^31 



(43) 



-E 3 i/T d 



1 + £42 



£31 



lO^cm" 3 we get AT„ 



130 K, substituting Nh 2 = 3 x 
, = 3.5 K. 



3.3. Comparison to the de Jong model 

Let us compare the de Jong model to the cold dust model. 
Using the escape probability approximation, the radiation 
field is now J = (1 — p)S, where p is the probability for 
photons to escape without interactions and we neglected 
here the influence by the dust. The source function is then 



S = 



e + p 



-B. 



(44) 




Fig. 5. Contour plots of the levels of the constant 6i6~523 
maser optical depth r m (solid curves) at the plane gas-dust 
temperature difference - gas-to-dust mass ratio for the follow- 
ing set of parameters Nh 2 = 10 9 cm~ 3 , T d = 130 K, / d = 0.01, 
and H = 10 14 cm. The dashed curves show the corresponding 
contours if the de Jong mechanism is not operating (i.e. the 
photon escape from the slab is neglected and the K2 and L2 
terms in equation (12) are omitted). 

In the optically thin case, p — ► 1, and the source function 
is then identical to that in the case of the large amount of 
the cold dust (B d = 0, see cq. [40]). Thus, in the optically 
thin limit the de Jong model is similar to the (very) cold 
dust - hot gas model. 

4. Results 

4.1. Maser transition 6i6~523 

The main results of our calculations for the case of amor- 
phous ice of the size of 0.1 fim are shown on Figs. 4-6. The 
dependences of the inversion and the maser optical depth 
on iVn 2 and fn 2 o/fd for the fixed dust and gas tempera- 
tures are shown in Fig. 4. One sees that when the temper- 
ature difference is small (Fig. 4b), the inversion appears 
in the region where -/VH 2 /H 2 o//d ~ 10 6 within two orders 
of magnitude. The inversion is weak for small /h 2 o//<i- 
We tested the influence of the de Jong mechanism on the 
results in this region by neglecting the terms K% and L2 in 
the escape probabilities (12). In that case (dashed curve 
in Fig. 4b), the inversion is absent there because a higher 
temperature difference is needed (see eq. [43]) to invert 
the populations. 

When AT is larger (see Fig. 4c, d), the inversion be- 
comes larger and the inversion region increases in size, 
spreading also over to the region with a small water con- 
tent (/h 2 o//c! < 10~ 6 ). Here the radiation field is com- 
pletely dominated by the dust and the inversion depends 
on the collisional rates defined by N^ 2 only (for the fixed 
T). At lower dust temperature T d = 30 K (Fig. 4a), the 



Fig. 6. Temperature dependence of the maser optical depth for 
the same six sets of parameters (Nu 2 , fn 2 o/fd) as in Fig. 2. 
We fixed H = 10 14 cm, T d = 130 K, f d = 0.01. The solid curve 
corresponds to analytical formula (45) with free normalization. 

optical depth becomes larger by about 20% comparing 
with T d = 130 K (Fig. 4c). 

We should point out that the inversion disappears at 
./Vh 2 ~ 10 12 cm -3 because of the level thermalization by 
collisions. This is lower than our analytical estimate (42) 
which neglected many collisional transitions. On the right 
side the inversion region is bounded by A^H 2 /H 2 o//d < 
10 8 ' 5 for T = 250 K, which is very close to our estimate 
(37). 

Let us note that the maximum inversion occurs when 
the water-to-dust mass ratio /h 2 o//(J ~ 10 5 /A r n 2 . The 
maser optical depth, on the other hand, is proportional 
to Nn 2 oAn m which reaches the maximum at fu 2 o/fd ~ 
10 7 ' 5 /W H2 for T = 250 K. This corresponds to 7V H2 o = 
10 6 ' 5 /d- The inversion disappears at Nn 2 o > 10 7-5 /d be- 
cause the amount of dust (comparing with water) is not 
sufficient here to produce the inversion. Thus for a given 
dust content there exists an optimal water concentration 
that produces the strongest maser. In order to calculate 
the maser intensity (sec cq. [23]) or the brightness temper- 
ature in the line center Tb r ~ |T ox | exp(r m / ^/tt), one needs 
to have a rough estimate for the excitation temperature. 
Because the dependence is linear, the error in that is not 
important. At the maximum of r m , the excitation tem- 
perature varies between —10 and —100 K when T varies 
between 150 and 500 K. 

On Fig. 5 we present the dependence of r m on fn 2 o/fd 
and AT for the fixed value of Nn 2 — 10 9 cm~ 3 . We see 
that the inversion disappears at large /H 2 o//d because 
the levels are then thermalized at the gas temperature. 
At small /H 2 o//d, the medium is optically thin for the 
line and the dust emission (rj < 1 for most infrared lines 
at this A% 2 for /<j = 0.01). Here we also tested the influ- 
ence of the photon escape from the surface. The dashed 
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Fig. 7. Contours of the constant maser optical depth for the strongest masing transitions of ortho- water for T<j = 130 K, T = 400 
K. Other parameters are the same as in Fig. 4. 



curves show the maser optical depth when the escape is 
neglected. These results show that at small AT the inver- 
sion is produced by the de Jong mechanism, while without 
it the minimum AT to invert the population is about <~ 1 
K. For small fn 2 o/ fd, however, the inversion appears only 
for AT > 15 K. At higher N H2 ~ 10 10 c m - 3 this limit is 
3.5 K (see dashed curve in Fig. 9) which is very close to 
the predictions of our four- level model (43). 

Fig. 6 shows the dependence of the optical depth in 
the main masing transition r m on the gas temperature 
for the same six sets of parameters (Nn 2 , /h 2 o//ci) &s in 
Fig. 2. One sees that in most points the inversion appears 
when gas temperature exceeds Tj. At point 4, however, 
the inversion exists even at smaller T. This results from 
the action of the de Jong mechanism in an optically thin 
(in line) medium for a small dust content (see also Fig. 3). 



The optical depth increases sharply at small AT and 
saturates at T <~ 500 K. This behavior is accurately de- 
scribed by our analytical formula (32). The optical depth 
is thus proportional to 

r m oc (c- E ^' T - e-WT d ) (45) 
where the VT factor comes from the Doppler width. 
4.2. Other maser transitions 

Our simulations show that in addition to the 6i6 — > 523 
maser, there appear many other masers. The strongest 
masers for ortho- water are listed in Table 1 and for para- 
water in Table 2. The levels of the constant optical depth 
are shown in Figs. 7 and 8, respectively. Some of these 
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Fig. 8. Same as in Fig. 7, but for para-water. 



masers were found in the calculations of Chandra ct al. 
(1984a) and Wallin & Watson (1997). Most of them ap- 
pear also in models based on the de Jong (1973) mecha- 
nism (see also Cooke & Elitzur 1985; Neufeld & Melnick 
1991). The 4 14 — > 3 2 i and 3i 3 — > 2 2 o transitions have been 
detected in the Orion cloud by Phillips et al. (1980) and 
Waters et al. (1980). The masing emission in IO29 — * 936 
line was observed in a wide variety of sources where 1.35 
cm maser was detected (Menten, Melnik, & Phillips 1990). 

We show the isolines of the zero inversion of the level 
populations as a function of the gas temperature for the 
fixed iVn 2 in Figs. 9 and 10. The inversion of the level pop- 
ulations appears on the right side from the corresponding 
curves. We see that the inversion in most of the transi- 
tions arises when the difference between the gas and dust 
temperature is larger than 30-100 K and one does not 
expect to see a significant signal in other water masing 



transitions if AT is small. This is related to the fact that 
the energy separation between the corresponding levels is 
much larger than that for the 616 — > 523 maser. Many 
masers also require the water concentration two orders of 
magnitude smaller than that where the 616 — > 52 3 and 
IO29 ~ > masers are strongest (see Figs. 7 and 8). In 
these respects the cold dust - hot gas model significantly 
differs from the de Jong model where many masing tran- 
sitions appear simultaneously. 

4.3. Other dust types and sizes 

Since the optical properties of the ice in the far infrared 
spectral range for different grain size 0.01 /tni < a < 1 /im 
are very similar, the maser strength is almost indepen- 
dent of a. The crystalline ice has a stronger peak in k<j at 
45 /im and therefore the 616 — 52 3 maser optical depth is 
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TABLE 1 
Ortho-water 
masing transitions 



# Transition A fim 



1 


3l2 


-> 2 2i 


260 


2 


4l4 


-> 3 2 1 


789 


3 


4 2 3 


— > 330 


669 


4 


6l6 


— > 5 23 


13475 


5 


532 


-> 4 4 i 


483 


6 


625 


— > 532 


227 


7 


634 


— > 541 


259 


8 


643 


— * 5so 


683 


9 


827 


^7 3 4 


231 


10 


743 


— > 652 


235 


11 


752 


— > 6 6 i 


677 


12 


1029 


-> 9 36 


933 



TABLE 2 
Para-water 
masing transitions 

# Transition A /Km 



1 


3l3 




2o2 


1636 


2 


5l5 




4 2 2 


922 


3 


524 




4 3 i 


309 


4 


533 




606 


632 


5 


633 




542 


194 


6 


726 




633 


208 


7 


642 




5si 


637 


8 


735 




642 


170 


9 


744 




651 


256 


10 


928 




835 


331 


11 


844 




753 


139 


12 


853 




762 


252 




1 10 100 

T - T d [K] 

Fig. 9. Isolines of the zero inversion of the strongest masing 
transitions for the same parameters as in Fig. 5. The labels on 
the contours correspond to the numbering in Table 1. Panel (a) 
is for Nu 2 = 10 10 cm -3 ; (b) corresponds to Nn 2 = 10 9 cm~ 3 . 
The inversion of the maser level populations exists on the right 
side of the corresponding isoline. The left dashed curves show 
the zero inversion for the 616 — > 523 maser when the de Jong 
mechanism is turned off. 



30% larger than that for the amorphous ice. We repeated 
calculations with 0.1 /Km silicate and graphite grains. The 
optical depth decreases by 15% and 25%, respectively. 

We also computed the maser conditions for the mix- 
ture of the graphite and SiC dust. We consider the size 
distribution of the dust grains dn(a) oc a~ 15 da, which 
extends from a m j n = 0.005 ^im to a max = 10 /xm. To al- 
low easy comparison with the previous results (Collison & 
Watson 1995; Wallin & Watson 1997), we use an approx- 
imation of the results by Laor & Draine (1993, see their 
Fig. 6) for the dust cross-section. We assume the dust 
opacity in the form «d = 600(/d/0.01) for A < 50 jum and 
K d = 600(/ d /0.01)(50 ^m/A) 2 [cm^" 1 ], for A > 50 /zm. 
The resulting maser optical depth is very similar to that 
for the crystalline ice. 



4.4. Approximating the inversion efficiency for 
616 — ► 523 maser 

The maser luminosity is a function of the optical depth r m 
which depends on the inversion An m . In many astrophys- 
ical problems, one would like to make simple estimations 
of the inversion and optical depth not repeating cumber- 
some calculations of the water molecular level populations. 
Thus one would like to have simple analytical approximate 
formulae for this purpose. 

It is possible to design a formula for the inversion that 
describes it with an error of a factor of two using the 
four-level model from Sect. 3. We can further simplify the 
model by assuming that the upper masing level 616 is pop- 
ulated directly from the 4i4 level. We then arrive to a sys- 
tem of two equations, each corresponding to the two-level 
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Fig. 10. Same as Fig. 9, but for para-water. The labels on the 
contours correspond to the numbering in Table 2. 

model. We now can introduce a pseudo-transition 1 — > 4, 
and prescribe to it some collisional and radiative rates. 

We propose the following formula to compute the in- 
version: 



An m = 



n-i ( S41 



S31 



(46) 



gi \ 1 + 1 + S31 

The population at the lower level 1 is given by the follow- 



ing expression 



91 



= 0.75C/- 1 rt exp(-289/T), 



(47) 



where, based on the results of simulations, we can approx- 
imate the partition function in the range 10 < T < 500 K 
as 



C/ part = 2.86[l + (T/42) 1 - 4 ]. 



(48) 



Approximating the dependence of the collisional rate on 
the gas temperature by a power-law, we get the ratio 

0.92 



T 



e 3 i - 0.033^ I — ) (1 - exp[-319/T]) . (49) 



The dust influence is parametrized by (3: 
/3 3 i = 2 x W-^T/WO^if^o/U)- 1 , 



(50) 



and £31 is computed using equation (14). We now can 
assume that the 4 — > 1 transition is characterized by 
£41 = l.le3i and /?4i = 0. 6/331. The source functions are 
found from equation (29). These formulae give a rather 
good agreement with the results of numerical simulations 
(within a factor of two) . 

5. Conclusions 

We consider the cold dust - hot gas mechanism of pump- 
ing of the water masers. To obtain the inversion of the 
maser level populations we take into account the first 45 
rotational levels of ortho- and para-H^O molecules and 
all possible collisional and radiative transitions between 
them. In the case of the large optical depth in the main 
pumping lines and the unsaturated masers, the radiative 
transfer problem can be significantly simplified, allowing 
us to investigate the maser efficiency for the set of the 
main parameters, such as the hydrogen concentration, the 
gas-to-dust mass ratio and the gas temperature. 

As it was suggested by Deguchi (1981), the pump- 
ing cycle 4i4 — ► 5os — ► 616 — > 523 - * 4i4 mainly deter- 
mines the inversion of the 616 — 523 maser levels popu- 
lations. We use this four-level model for the analysis of 
our numerical results. We also suggest approximate for- 
mulae for the dependence of the inversion on the hy- 
drogen concentration, the gas-to-dust mass ratio and the 
gas and dust temperatures, that could be used for mod- 
eling the astrophysical sources. We find that the inver- 
sion is the largest when the water-to-dust mass ratio is 
/H 2 o//d ~ 10 5 /Ah 2 , an d the maser optical depth is the 
largest when Ah 2 o ~ (10 6 + 10 7 )/ d - 

The inversion of the maser level populations as a func- 
tion of the same parameters is calculated in the case of the 
amorphous and crystalline ice, silicate as well as graphite 
dust. The results of these calculations are very similar. 
Because the inversion depends on the dust mass density, 
our results can be applied to any size distribution of the 
dust particles. 

Our analysis shows that in the optically thick environ- 
ments the inversion in the 616 — > 523 transition appears 
when the gas is just ~ 1 K hotter than the dust, while 
most of other transitions start to be inverted at much 
larger AT ~ 30—100 K. Most of them require also smaller 
water-to-dust mass ratio /h 2 o//c1- This is very different 
from the predictions of the de Jong model where many 
transitions are inverted simultaneously. If the cold dust - 
hot gas model is the correct one, the relative ratios of the 
luminosities in different masing transitions could provide 
constraints on the gas and dust temperatures, water-to- 
dust mass ratio, and hydrogen concentration. 
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